Non-adiabatic many-atom quantum state control in few-well systems 
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We present a fast scheme for arbitrary unitary control of interacting bosonic atoms in a double- 
well. Assuming fixed inter-well tunnelling rate and intra-well interaction strength, we control the 
many-atom state by a discrete sequence of shifts of the single-well energies. For strong interactions, 
resonant tunnelling transitions implement beam-splitter U(2) rotations among atom number eigen- 
states, which can be combined and, thus, permit full controllability. By numerically optimizing such 
sequences of couplings at avoided level crossings (CALC), we extend the realm of full controllability 
to a wide range of realistic interaction parameters, while we remain in the simple control space. We 
demonstrate the efficiency and the high achievable fidelity of our proposal with non-adiabatic pop- 
ulation transfer, iVOOiV-state creation, a C-NOT gate, and a transistor- like, conditional evolution 
of several atoms. 

PACS numbers: 05.30.Jp, 03.75.Lm, 05.60.Gg, 37.10.Jk 



The possibility to count and manipulate individual 
particles makes ultracold atoms in optical lattices a 
highly attractive system for quantum state engineering 
and quantum computation [1-4 . While the fundamen- 
tal building blocks for the manipulation and detection of 
individual atoms have been established, the control over 
more complex quantum states in ultracold-atom systems 
remains a challenge [5]. 

In the paradigmatic double-well system, proposals ex- 
ist for certain tasks such as the full transfer of atoms be- 
tween the wells [5] , and the creation of particular classes 
of entangled states [5J [7] , TVOOiV-states [5] and squeezed 
states [§]■ Arbitrary state control of the atomic occu- 
pation dynamics is a desirable key component to realize 
quantum atomtronics |10H17j . but no universal protocol 
has so far been proposed for this endeavour. 

General proofs of controllability have been achieved 
from two different perspectives: Any iV-dimcnsional 
quantum system can always be controlled by the de- 
composition of the desired unitary U(iV) operation in 
elementary U(2) blocks 18-20 . On the other hand, 
every kinematically allowed transformation can be re- 
alized via "bang-bang" control [3TJ, i.e. by alternat- 
ing the time-evolution induced by two sufficiently non- 
commuting Hamiltonians |22j . The existence of a con- 
trol scheme does not, however, provide a solution to the 
practical task of reliably and quickly finding simple and 
robust control sequences. 

In this letter we achieve perfect unitary control by com- 
bining the two above approaches. We use a numerically 
optimized sequence of couplings at avoided level crossings 
(CALC) to bridge the gap between abstract existence 
proofs for control sequences [3TJ [55] and their actual, re- 
liable and robust implementation in many-particle few- 
well systems [?H9"ll23|. In particular, we focus on the truly 
high-fidelity (> 0.999) preparation of quantum many- 
body states as possible resources for quantum informa- 
tion technologies. As examples, we apply our approach 



to the non-adiabatic transport of N atoms from one well 
to another, to the creation of A^OOA^-states, and to the 
controlled dynamics of a single-species atom-transistor. 
Randomly sampling over target states reveals that truly 
every quantum state can be created within a wide range 
of realistic parameters. 
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FIG. 1: (color online) (a) Bosonic atoms in a double-well sys- 
tem with adjustable tilt e. (b) Fock-state energy levels as a 
function of the tilt e in the limit U 3> J, with highlighted first- 
order (red solid) and higher-order (black dotted) avoided level 
crossings, (c) Two-state couplers at avoided level crossings 
(large orange arrows) and phase-shifters at well-separated en- 
ergies (blue dotted lines), (d) Time-dependent tilt e(i) that 
realizes the couplings and phase shifts in the 4-dimensional 
space of N = 3 particles, as shown in (c). 



We consider a double-well system with N bosons, il- 
lustrated in Fig. [lja), which is described by the Bose- 
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Hubbard Hamiltonian, 

#bh(c) = -~(a\a 2 + a\a-^ (1) 
t* ,\ U i- 

3 = 1 

where U is the collisional interaction strength, J is the 
inter-well tunnelling strength, and e is the double-well 
tilt, proportional to the difference between single-particle 
energies in the wells. We assume that e can be varied in 
the experiment [24J [25]. In the case of U > J the N + 1 
two-mode Fock-states 

\N,0),\N-l,l),...,\Q,N), (2) 

constitute eigenstates of Hbh(^) for most values of e (see 

Expanding Eq. (JlJ in the Fock-state basis, one read- 
ily finds that, for U ^ 0, two Hamiltonians associated 
with different values of e fulfill the requirements of "bang- 
bang" -control [23112]; such that alternating applications 
of the two Hamiltonians can provide any unitary oper- 
ation. Already for N = 3, however, the restriction to 
only two values of e is unnatural and inefficient, and in 
the following we shall suggest a simple control sequence 
that exploits the possibility for e to attain a wide range 
of values. 

We start by giving a systematic way to obtain the 
CALC sequence for strong interactions, U S> J, for which 
the energy levels are depicted in Fig.[]|b). This sequence 
will then be used as a starting point for numerical opti- 
mization for realistic values of U / J. Two states \n, N—n) 
and \n + m, N — n — to) are degenerate when 

e = Vn,m ■= 2 «- N + m (n = 0, . . . , N - 1). (3) 

Correlated tunnelling of to > 1 particles can only oc- 
cur off-resonantly and it is thus negligible in the regime 
U ^> J. Instead, single-particle tunnelling (m = 1) dom- 
inates the picture: Given the tilt e = r/ n .i, the Hamil- 
tonian [Eq. ([I])] couples a degenerate pair of states with 
the bosonically enhanced frequency 

ui n = J\Jn + ly/N - n, (4) 

while possible accidental degeneracies with m > 1 can be 
neglected and the Hamiltonian remains diagonal in all 
the other Fock-states. This immediately suggests that 
all atoms can be transferred from the right to the left 
well when the tilt is brought quickly from large negative 
values to 770.1, then 7711 etc., spending the time 

i„ =ir/uj ni (5) 

at each value T) n 1 to ensure complete state transfer at the 
avoided level crossing. In the limit N — ¥ 00, we recover 
the continuous transfer protocol derived in Ref. [5J. 



Staying in the regime U ^> </, a given duration spent 
at a first-order avoided level crossing leads to the equiv- 
alent of a beam splitter operation with a certain reflec- 
tivity on a pair of optical modes. A far-reaching result 
in quantum optics is that a general TJ(N + I) multi- 
port beam splitter can be implemented by an arrange- 
ment of normal N(N + l)/2 two-mode beam splitters 
and (N + 1)(N + 2)/2 phase shifters [19]. Here, we re- 
quire a maximum of N — n operations on each pair of 
states \n, N — n) and \n + 1, N — n— 1), following the ex- 
plicit decomposition of Ref. [T!5] ■ Exploiting these results 
for our system, we find that a CALC sequence that is 
constructed by strictly following the beam splitter anal- 
ogy can implement every unitary operation in the limit 
U ^> J. Assuming that phase shifts, acquired at a rate 
proportional to eU, take a vanishing fraction of the time, 
we can give a bound on the required total protocol dura- 
tion T via 

The validity of the above beam splitter analogy heavily 
relies on the assumption U J. This physical regime is, 
however, not desirable in the laboratory, since inelastic 
collisional losses become important for strong interaction 
U, whereas the use of very small tunnelling strengths J 
prohibitively increases the overall time-scale [Eq. ([6|] of 
the process. At realistic, finite values of U/J ~ 5 ... 10, 
the avoided level crossings in Fig. [ljb) become broad, 
such that two-mode Fock-states cease to be eigenstates 
of the problem. Additionally, higher-order avoided level 
crossings also become relevant. Consequently, by setting 
e = rj n ,1, we not only induce transfer between \n,N — n) 
and |n+l, N—n—1) but also between other pairs of states 
that differ by one or several particles. The simple picture 
of CALC dynamics as a sequence of beam splitters and 
phase-shifters is therefore not valid anymore. 

For finite U/J, an adiabatic change of the tilt from 
large negative to positive values induces the transition 
[0, N) — > |7V, 0) , via the slowly changing lowest energy 
eigenstate of the system. The total protocol duration, 
however, is then much larger than 1/J. In the follow- 
ing, we shall numerically optimize CALC sequences to 
account for finite U / J and obtain a high-fidelity scheme 
that is universal and practical, in the sense that any tar- 
get state I ^target) ca n be reached starting with a given 
initial state | ^initial) in few, simple operations. 

A general control sequence consists of M values of the 
tilt ej (1 < j < M), which are each applied for a duration 
tj . The sequence implements the unitary evolution 

M 

U[tA = J]V^BH( £fc )t fc; (J) 

fc=l 

where the product is understood to respect the time- 
ordering between the pulses. The fidelity of the state 
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FIG. 2: (color online) (a) Target-state infidelity 1 — .F for full 
population transfer |0,JV) — > \N,0) without optimization, as 
a function of the interaction strength U/J, for TV = 2, 3, 5, 10 
(red dotted, green dot-dashed, black solid, and blue dashed 
line, respectively). Optimizing numerically always yields ar- 
bitrarily high fidelity, for any value of U/J (not shown), (b) 
Total protocol duration for the numerically optimized popu- 
lation transfer. The thin dashed lines correspond to the ideal 
total protocol duration in the limit U/J—¥ oo, T ■ 
with tk given by Eq. (jSJ. 



preparation is quantified by the overlap of the prepared 
state with the desired target state, 



T= | (*targot|W[^e| |*i„itial} 



(8) 



The optimal vectors of time and tilt variables (t, e) 
will not exactly coincide with the ones identified for 
U/J — > oo using the decomposition given in [12], but we 
find it convenient to use the CALC sequence as an ini- 
tial guess for numerical optimization in which the pulse 
lengths tj as well as the tilts ej are adjusted, while the to- 
tal number of steps M is kept constant. For finite U/J, 
perfect fidelity can be achieved in most cases via such 
numerically optimized CALC sequences. When no satis- 
factory solution is found, the number of steps M can be 
increased to extend perfect fidelity to the regime of even 
smaller U /J. 

As a first application, we discuss the full transfer of N 
particles between the wells, for which the unoptimized 
CALC sequence with M = N resonant tilt values that 
are applied for durations given by Eq. Q only leads to 
acceptable fidelities when the interaction U/J is large 
and the number of particles small, as shown in Fig. [2ja). 
Optimizing the times t k , or, alternatively, the tilts e k 
does not lead to a satisfactory solution, but by simultane- 
ously optimizing both sets of parameters using a standard 
Nelder-Mead simplex algorithm |22], we quickly reach, 
for any value of U / J, an infidelity that is only constrained 
by numerical accuracy. We show the evolution of the 
wave-function components for N = 3 in Fig. [3] In (a) , 
the interaction U/J = 40 leads to a high fidelity, with- 
out optimization of the CALC sequence. Using the same 
sequence, however, an unacceptable, low fidelity is ob- 
tained for small interactions, U/J — 1, in panel (b), but 
by optimizing, high fidelity is again recovered combin- 
ing single-particle and correlated tunnelling through dif- 
ferent intermediate states [panel (c)]. The unoptimized 




FIG. 3: (color online) Optimized and non-optimized se- 
quences and population of state vectors for the full transfer 
of TV = 3 particles from the left to the right well, starting 
with 1 3, 0} . (a-c) Population of the different Fock-states and 
(d) pulse sequence, as a function of time, (a) The interaction 
U / J = 40 is large, such that a non-optimized control sequence 
leads to a high fidelity, (b) The same control sequence, how- 
ever, transfers the atoms only with around 68% fidelity when 
the interaction is small, U/J — 1. (c) By optimizing the 
control sequence numerically, perfect fidelity can be achieved 
again, (d) The optimized pulse-sequence (dashed blue) differs 
from the analytically obtained one (solid red). The total time 
required for the optimized transfer is only slightly larger than 
for the non-optimized sequence. 



(solid red) and optimized (dashed blue) CALC sequences 
[panel (d)] are very similar, and as shown in Fig. [2^b), 
the total protocol durations of the unoptimized and of the 
optimized sequences do not deviate considerably. The fi- 
delity of the process depends on the accuracy of the time 
and tilt settings; F > 0.999 is reached as long as the 
times tk and the tilts e k are subjected to an error up to 
1% and 0.01, respectively. 

Although, for U = 0, the trivial solution to the trans- 
port problem consists in setting eo — for a time 
to = ir/ J so that all particles tunnel at the same time, 
this protocol is not robust against small deviations from 
[7 = 0, and even for small values of U/J, the CALC se- 
quence provides a better starting point for optimization 
than the singular solution at U = 0. 

While unconstrained numerically optimized CALC se- 
quences are efficient for any U/J, the "bang-bang" ap- 
proach with only two different values of e leads to con- 
siderably more complicated protocols: For instance, we 
needed 40 steps to achieve J- > 0.99 for the full transfer 
of N = 5 particles at U/J = 5, and the number of steps 
increases further for larger U/J. 

As an example of a non-trivial transformation, we con- 
sider the complete transfer protocol for N atoms, but we 
maintain the first tilt e = 770,1 for half the transfer time, 
to = tt/(2lu ), while the other tilts are kept as in Eq. (5). 
As can be readily understood from Fig. [ljb), this se- 
quence of N steps transforms the Fock-state |Q, TV) into 
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FIG. 4: (color online) Target-state infidelities 1 — T for 
/VOOiV-state creation, |0,iV> -> 1/v^ (|JV, 0) + e^|0, TV)) , 
without and with optimization, as a function of the interac- 
tion strength U/J. Blue dashed line: No optimization. Solid 
black line: Numerically optimized, for M — N. Red dot- 
ted line: Numerically optimized, for M = IN. To facilitate 
reading, we show the worst fidelity obtained for any interac- 
tion smaller than U/J after optimization, which renders the 
curve monotonic. The achievable infidelity after optimization 
abruptly changes from an arbitrarily low value (here: 10™ 8 , 
not shown) to a finite value at a certain U/J, depending on 
the number of particles. 



the entangled AOOA-state, (\N,0) + exp(^6)|0, N))/y/2, 
where the relative phase cj> can be adjusted by the time 
spent at tilt values away from the avoided level cross- 
ings. In Fig. |4j we display the results for the generation 
of AOOA-states. The dashed blue line represents the in- 
fidelity for the unoptimized CALC sequence, the solid 
black line shows the result after numerical optimization. 
In contrast to the complete transfer of A atoms between 
the wells, the infidelity of the optimized CALC sequence 
with M — N steps rises abruptly below a critical value 
of U/J. This is consistent with the geometrical picture 
of SU(2) rotations around non-orthogonal axes [27| 128]. 
and a mere doubling of the number of time-steps allowed 
in our optimization, M — > 2M, significantly extends the 
range of U/J that allows high fidelity AOOA-state cre- 
ation (dotted red line). 

In general, we find that every two-mode state of A 
atoms can be generated by numerically optimized CALC. 
This universality was confirmed by picking a large ran- 
dom set of states that were all successfully generated for 
finite U/J. 

Finally, we apply numerically optimized CALC to a 
three-well single-atom transistor: Conditioned on the 
presence of an atom in the middle well of a three-well 
potential, all atoms in the left well are transferred to 
the right well. We assume that the single-particle energy 
in the right (left) well can be shifted with respect to the 
other two wells in the first (second) step of the procedure. 
The process can be implemented by two subsequent com- 
plete population transfer sequences that are tailored such 
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FIG. 5: (color online) Transistor-like time-evolutions, (a) 
Without optimization, the transistor sequence for TV = 
5, U/J = 4 only yields T ~ 0.73: The |5, 0)-component (black 
dashed) is not eventually populated with unit probability, and 
a finite population stays in the |5, f)-state (blue solid), at the 
expense of the |0, 6}-component (red solid). Other compo- 
nents are drawn as thin gray (N + f-particle sequence) and 
black dashed lines (TV-particle sequence), (b) By optimizing 
the control sequence with 40 steps (M = 8N), T > 0.999 
is reached, (c) Infidelity of transistor-sequence for different 
numbers of particles and U/J — 4, as a function of the num- 
ber of iterations applied; denotes no optimization, otherwise 
the number denotes the relative increase in the number of 
steps, M/N. (d) Non-optimized solution for U ^> J (dashed 
blue) and optimized pulse-sequence for U / J = 4 (solid dark 
red) . A short spike at t ~ extending from e = —45 to e = 38 
was omitted in the figure. 
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| A, 1,0) 
| A, 0,0) 



|0,A + 1,0) 
| A, 0,0) 



(») 



|0,1,A), 
| AT, 0,0). 



(9) 



The second step (ii) is a special case of our double well 
control, and it only needs to be implemented successfully 
for the A + 1 atoms present in the middle well in the first 
line, whereas the first step (i) needs to be optimized to 
ensure that transfer only occurs into the already occupied 
middle well. We assess the control sequence by the prod- 
uct of the fidelities of the two processes [Eq. (|9|]. Due to 
the complexity of the dynamics, optimizing the CALC se- 
quence while keeping the number of control steps M = A 
does not yield acceptable fidelities. By iteratively dou- 
bling the number of steps M while keeping the total pro- 
tocol duration approximately constant, we do, however, 
reach arbitrarily high fidelities for finite U / J. For A = 1, 
the evolution [Eq. ^] can be used as the main building 
block of a two-qubit C-NOT-gate. We reach T > 0.999 
for all U/J > 1 with only M < 8 settings of the tilt 
variables. For A = 5 and U/J — 4, we show the unopti- 
mized and the fully optimized time-evolution of the rele- 
vant components of the wave- function in Fig.[5]Ja,b), and 
the evolution of the achievable fidelities with the num- 
ber of iterations in Fig. [5jc). The minimal interaction 
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strength required for T > 0.999 at a constant number of 
steps M increases with the number of particles, just like 
for the 7V0(W-state creation above. Due to the complex- 
ity of the dynamics, errors in the time and tilt settings 
jeopardise the process fidelity more severely than for the 
two-well dynamics; T > 0.99 is obtained for times (tilts) 
that are accurate at the level of 1% (0.005). We stress 
that unlike other atomic transistor proposals [101 Hi] our 
approach does not involve two different atomic species. 

Although any system with sufficient coupling between 
states permits full control |HJ[52], it is in many situations 
not clear how to construct efficient and simple control se- 
quences. In few-well systems with interacting bosons, the 
use of CALC sequences is a natural strategy that exploits 
the analogy to beam splitters in the limit of strong in- 
teractions. Numerically optimized CALC sequences can 
then achieve complex conditional dynamics, as required 
for quantum computation and atomtronics. The CALC 
sequence for U/J 3> 1 can be easily obtained [19], it typ- 
ically provides a fidelity above 90% for moderate values 
of U, which, on the one hand, is insufficient for appli- 
cations, but, on the other hand, provides a good start- 
ing point for numerical optimization to reach a fidelity 
that is arbitrarily close to unity. Our choice of optimiza- 
tion parameters is the simplest extension of the original 
CALC sequence. Thus, thanks to the manageable, re- 
stricted number of optimization parameters, control via 
numerically optimized CALC is also feasible for complex 
Hamiltonians 29] . 

A natural extension is a protocol that tolerates particle 
number fluctuations: As long as the constraints of unitar- 
ity are not violated, there is no natural boundary for the 
design of complex conditional control sequences. Simi- 
larly, sequences that tolerate larger errors in the time and 
tilt sequences are desirable. From a control-theory per- 
spective, an extension of the SU(2) theory of generalized 
Euler-angles for non-orthogonal rotation axes [27J HH] to 
larger spaces may complement our numerical results and 
give analytic insight in the minimal interaction strength 
U/J that is required for perfect control under the con- 
straint of a limited number of steps M. An interesting 
extension of the scheme described by Eq. Q would be the 
adaptation to more sophisticated methods for "shortcuts 
to adiabadicity" |30) . 
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